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NOTE  ADDED  DECEMBER  1976 


This  work  was  done  early  in  1975,  and  this  report  was  written  in  the 
summer  of  that  year.  Because  of  the  general  sensitivity  of  the  subject,  the 
report  received  only  very  limited  distribution  to  interested  ERDA  representatives 
at  that  time.  In  fact,  we  were  asked  to  withhold  further  distribution  until  addi- 
tional calculations  were  completed  and  a more  mature  understanding  of  the  in- 
fluence of  material  strength  on  residual  stress  in  layered  media  was  obtained. 

Most  of  the  requested  additional  calculations  have  been  finished,  and 
many  of  the  results  have  been  presented  to  an  Earth  Motion  Calculator's  meeting 
held  in  La  Jolla  on  May  4,  1976.  Since  the  multidimensional  calculations  re- 
quested at  that  time  have  now  been  done  and  reported  by  LLL  at  the  66th  CEP 
meeting,  it  seems  appropriate,  for  archival  reasons  at  least,  to  distribute  the 
earlier  report. 

One  minor  change  has  been  made  in  this  document;  namely,  Reference  12 
has  been  changed  to  a more  complete  and  up-to-date  presentation  of  the  influence 
of  material  properties  on  residual  stress  and  cavity  radius.  As  mentioned  above, 
additional  and  better  calculations  of  the  sort  reported  here  have  been  completed, 
and  a few  more  are  to  be  done.  An  extended  version  of  this  report  will  be  distri- 
buted in  due  course. 


ttCBSIONfor 


NTtS  Witte  Section 

DOC  Butt  Section  □ 

UNANNOUNCED  □ 

JUSTIFICATION 


vf 


BY  - 

BisnuEinir^/.^vM1  ASiitiv  coses 

Dist.  ••  j(/c>f  SPECIAL 


A 


■ "V*  b- 


7;7v/\v 

- ' ....  - i>\  •-  • - 


1.  INTRODUCTION 


This  report  discusses  a recent  study  at  Systems,  Science  and  Software 
of  Baneberry,  an  underground  nuclear  event  detonated  in  Yucca  Flat  in  December 
1970.  Baneberry  failed  to  contain  and  therefore  has  been  studied  in  great  detail. 
However,  a conclusive  answer  to  the  question  of  why  Baneberry  vented  is  still  not 
available. 

The  Baneberry  working  point  was  located  in  an  alluvial  soil  containing 
large  concentrations  of  montmorillonite  clay.  A layer  of  alluvial  material  having 
high  sonic  velocities  was  directly  above  the  working  point  while  a paleozoic  escarp- 
ment was  located  below  and  on  one  side  of  the  working  point.  Postshot  calculations 
by  Crowley  and  Cherry  (1971)  looked  into  the  effect  of  the  high  sonic  velocity  layer 
while  Cherry,  et  al.  (1974),  discussed  the  effect  of  the  escarpment.  The  general 
consensus  of  these  studies,  as  I understand  them,  is  that  the  containment  failure 
of  Baneberry  is  attributed  to  the  high  water  content  at  the  working  point  due  to 
saturated  clay,  which  resulted  in  extremely  high  coupling,  and  to  the  large  imped- 
ance contrasts  near  the  working  point,  which  resulted  in  enhanced  tension  failure 
caused  by  the  nonspherical  ground  motion  compared  with  other  events  in  Yucca  Flat. 

Since  Baneberry  remains  of  great  interest  today,  a near  field  ground 
motion  calculation  of  the  event  was  made  using  the  one-dimensional  Lagrangian 
code,  SKIPPER,  incorporating  the  most  recent  geological  and  material  properties 
data  available.  The  results  were  compared  with  similar  calculations  for  a Yucca 
Flat  dry  tuff,  described  in  Bache,  et  al.  (1974),  and  for  a Rainier  Mesa  saturated 
tuff  event  (Husky  Ace),  described  in  Rimer,  et  al  (1974).  A scenario  is  presented 
which  indicates  that  the  containment  failure  of  Baneberry  may  be  due  to  an  absence 
of  significant  compressive  residual  stresses  around  the  Baneberry  cavity.  Re- 
sidual stresses,  due  to  nonuniform  plastic  loading  of  the  rock,  are  observed  in 
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most  near  field  calculations  at  S l see  Rimer  (1976)]  and  reach  their  final  mag- 
nitudes at  approximately  the  time  the  cavity  attains  its  full  growth.  These  re- 
sidual stresses  make  the  cavity  produced  by  a nuclear  explosion  a relatively 
strong  structure  which  contains  gases  very  well.  In  the  absence  of  this  so-called 
"mystical  magical  membrane, " tensile  hydrofracturing  is  believed  to  have  taken 
place,  leading  to  venting  of  the  cavity  gases. 

The  plan  of  this  report  is  to  detail  in  Section  2 the  geological  and  ma- 
terial properties  data  and  modeling  which  were  used  in  the  calculations  mentioned 
above.  Section  3 discusses  the  results  of  these  calculations.  The  conclusions  of 
this  report  are  summarized  in  Section  4. 

2.  MATERIAL  PROPERTIES 

The  Baneberry  device  was  emplaced  at  a depth  of  912  ft  in  hole  U8d  on 
the  far  west  side  of  Yucca  Flat.  A 1971  postshot  analysis  described  the  working 
point  as  fully  saturated  with  an  abnormally  high  water  content  for  a WP  located  at 
least  900  ft  above  the  water  table  due  to  the  presence  of  large  amounts  of  clay 
(montmorillonite)  at  this  depth.  This  would  be  expected  to  give  greater  than  normal 
ground  shock  motion  when  compared  to  other  events  in  Yucca  Flat.  The  seismic 
yield  of  Baneberry  has  in  fact  been  estimated  postshot  (1971)  at  about  three  times 
the  device  yield  of  the  event,  while  the  measured  surface  motion  was  far  greater 
than  anticipated.  The  rode  above  the  WP  failed  to  contain  the  cavity  gases  which 
began  to  vent  to  the  surface  3.5  minutes  after  detonation. 

One-dimensional  calculations  were  made  postshot  by  Crowley  and 
Cherry  (1971)  to  study  the  effect  of  a high  sonic  velocity  layer  of  saturated  rock 
above  the  working  point.  Two  calculations  were  made,  one  for  a normal  site, 
modeled  as  highly  porous,  dry  alluvium,  and  another  in  which  230  ft  of  the  allu- 
vium directly  around  the  WP  was  replaced  by  the  high  velocity  layer.  The 


calculations  indicated  that  the  high  velocity  layer  results  in  the  same  surface  spall 
velocity  for  Baneberry  yield  as  a blast  of  five  times  the  yield  in  a normal  site. 


The  one-dimensional  calculation  of  Baneberry  presented  in  this  re- 
port utilizes  geologic  and  material  properties  data  not  available  in  1971  and  also 
continues  die  calculation  to  later  times  (a  few  seconds)  in  order  to  examine  the  re- 
sidual stresses  around  the  cavity. 


Preshot  material  properties  data  were  based  largely  on  a north-south 
line  of  satellite  holes  drilled  for  Discuss  Thrower  event  emplaced  in  U8a,  1500  ft 
north  of  U8d  augmented  with  some  cuttings  samples  from  U8d  itself.  The  material 
around  the  Baneberry  WP  was  determined  to  be  undisturbed  Ammonia  Tanks  tuff. 
Stephens,  et  al.  (1971),  have  studied  an  altered  tuff  from  the  1 114— ft  level  of  hole 
Ue8f,  approximately  600  ft  east  of  U8d.  From  a limited  amount  of  sample,  de- 
scribed as  Ammonia  Tanks  tuff,  they  obtained  both  pressure-volume  data  to  30  kbar 
and  a failure  surface. 

A new  exploratory  hole,  U8ei,  360  ft  south  of  U8d,  was  drilled  in 
September  1973  in  order  to  obtain  a better  description  of  the  geology  and  material 
properties  appropriate  to  Baneberry.  The  U.  S.  Geological  Survey  has  restudied 
the  cuttings  samples  in  U8d  and,  because  of  the  better  quality  of  samples  in  Ue8i, 
suggests  that  the  WP  may  be  in  colluvium  rather  than  Ammonia  Tanks  tuff.  Based 
on  Ue8i  samples,  Ramspott  (1975)  has  described  the  working  point  as  being  in  a 
layer  of  over  50%  montmorillonite  rock.  The  water  content  in  this  saturated  layer 
may  be  as  much  as  25%  by  weight.  A higher  overburden  density  than  usual  for 
Yucca  Flat  was  reported.  This  would  lead  to  a higher  initial  postshot  cavity  pres- 
sure. Ramspott  has  also  mapped  the  approximate  clay  content  of  the  rock  layers 
about  the  working  point.  In  the  immediate  vicinity,  there  is  some  justification 
based  on  clay  content  for  modeling  these  layers  as  spherical. 
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From  the  above  information,  which  is  often  conflicting,  a consistent 
set  of  material  properties  data  was  chosen  as  initial  conditions  for  a one-dimensional 
calculation  of  Baneberry  using  the  Lagrangian  SKIPPER  code.  The  initial  conditions 
consisted  of  a one-zone  spherical  cavity  source  containing  the  device  energy  at  zero 
time  and  three  spherical  layers  representing,  respectively,  the  high  clay  content 
working  point  layer  (50  ft  thick),  the  high  source  velocity  layer  (120  ft  thick),  and 
an  alluvial  layer  extending  to  the  surface.  This  modeling  should  be  considered  ap- 
plicable only  to  a vertical  section  through  the  working  point  to  the  surface.  No  at- 
tempt is  made  to  treat  the  paleozoic  escarpment.  With  the  exception  of  the  cavity 
source,  all  layers  are  modeled  using  the  Tillotson  equation  of  state,  the  p-a 
porous  crushup  model,  and  the  parabolic  failure  surface  described  in  the  Appendix 
of  this  report.  Table  A-l  of  the  Appendix  lists  the  material  properties  data  used 
in  the  calculation. 

For  the  working  point  layer  (material  1),  the  pressure-volume  curve 
was  based  on  P-V  data  from  Stephens,  et  al.  (1971),  while  the  yield  strength  was 
considered  negligible  due  to  the  clay  content.  This  is  believed  justified  since  the 
emplacement  hole  repeatedly  caved  in  during  drilling.  This  material  was  con- 
sidered fully  saturated.  The  alluvium  layer  (material  3)  was  based  on  equation 
of  state  data  for  Sedan  alluvium  obtained  from  C.  Hastings.  Air-filled  porosity 
was  computed  to  force  these  data  to  correspond  to  measured  density  for  Baneberry 
alluvium.  The  high  sound  speed  layer  (material  2)  was  judged  to  have  one-half  the 
yield  strength  of  the  alluvium  and  to  be  more  saturated  due  to  its  clay  content 
(greater  than  20%).  An  overburden  pressure  of  60  bars  was  included  in  the  modeling. 

The  cavity  source  region  was  chosen  to  be  initially  large  enough  to  just 
melt  70  metric  tons  of  rock  for  each  kiloton  of  device  yield.  The  cavity  pressure 
is  given  by  an  ideal  gas  law  where  gamma  is  a function  of  specific  volume 
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(Y=  1.03  + 0,9/  Pv  ).  This  expression  for  y was  obtained  for  tuff  modeled  as 
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silicon  dioxide  from  the  Oracle  Chemical  Equilibrium  Code  adapted  at  S from  the 
Tiger  Code  discussed  by  Cowperthwaite  and  Zwisler  (1971).  This  code  uses  the 
JCZ3  gas  equation  of  state  and  the  S3S  solid  equation  of  state  of  D.  Laird  (1975). 

The  results  of  this  Baneberry  calculation  will  be  compared  to  results 
from  calculations  for  an  average  Yucca  Flat  dry  tuff  described  in  Bache.  et  al. 
(1974),  and  for  a saturated  Rainier  Mesa  tuff  (Husky  Ace)  described  in  Rimer,  et 
al  (1974).  Since  material  properties  information  for  Yucca  Flat  events  is,  in  gen- 
eral, not  extensive,  an  average  dry  tuff  was  modeled  for  the  former  calculation, 
based  on  L.  Germain’s  (1974)  summary  of  data  for  all  dry  tuffs.  Using  average 
values  for  bulk  density,  grain  density,  water  content,  porosity,  and  saturation, 
an  equation  of  state  was  constructed  using  the  Tabular  Array  of  Mixture  Equation 
of  State  (TAMEOS)  scheme  described  in  Riney,  et  al  (1972).  This  scheme  mixes 
grain  density  rock  with  water  assuming  pressure  equilibrium  between  the  mixture 
and  its  components  in  order  to  generate  a Hugoniot  and  isent ropes.  The  resulting 
pressures,  energies,  and  specific  volumes  are  rearranged  into  a fast  table-lookup 
format.  A crushup  pressure  and  an  elastic  pressure  limit, needed  for  the  P-a 
porous  crushup  model,  were  calculated  from  empirical  formulas  available  in 
Butkovitch  (1973). 

Sufficient  material  properties  data  are  available  to  accurately  model 
the  almost  fully  saturated  tuff  found  at  Rainier  Mesa.  This  calculation  also  utilized 
the  TAMEOS  scheme  and  the  P-a  model.  A summary  of  the  data  used  in  the  three 
calculations  may  be  found  in  Table  A-l  of  the  Appendix. 

3.  NUMERICAL  RESULTS 

Results  from  a spherical  one-dimensional  Baneberry  calculation  were 
compared  with  earlier  results  from  calculations  for  a Rainier  Mesa  saturated  tuff 


and  a Yucca  Flat  dry  tuff.  Peak  velocity  versus  range  has  been  plotted  in  Fig.  1 
as  a measure  of  the  ground  shock  coupling  of  the  three  events.  Note  that  the 
Baneberry  free  surface  is  not  included  in  the  calculation.  All  curves  are  scaled 
to  the  Baneberry  yield.  Over  almost  the  entire  range  plotted,  Baneberry  peak 
velocity  is  intermediate  between  Rainier  Mesa  and  Yucca  Flat.  However,  Baneberry 
is  seen  to  couple  even  higher  than  Rainier  Mesa  in  the  first  layer  (r  < 1524  cm) 
containing  large  amounts  of  saturated  clay.  The  agreement  (shown  in  Table  1) 
between  the  Baneberry  calculation  and  experimental  measurements  obtained  from 
Crowley  and  Cherry  (1971)  is  excellent. 


Table  1.  Comparison  between  Baneberry  Data  and  Calculations 


Station  Location  (meters) 

Peak  Velocity  (meters/sec) 

Data 

Calculation 

130 

12.63 

12.5 

135 

10.47 

11.0 

Peak  radial  stress  vs  range  is  plotted  in  Fig.  2 for  the  same  three 
calculations.  As  in  Fig.  1,  Rainier  Mesa  saturated  tuff  plots  higher  than  Baneberry 
except  in  the  first  layer  of  highly  saturated  clay.  Baneberry  couples  much  higher 
in  peak  stress  than  Yucca  Flat  dry  tuff  in  the  first  two  layers.  However,  in  the  al- 
luvium layer  (r  * 5180  cm),  the  Baneberry  calculation  shows  lower  peak  radial 
stresses  than  the  dry  tuff  at  the  same  range.  This  result  follows  from  the  yield 
strengths  chosen  for  the  two  materials. 


Peak  Velocity  (cm/sec) 


Position  (cm) 


Figure  1 Peak  velocity  versus  range  for  the  Baneberry,  Rainier 
Mesa  and  Yucca  Flat  calculations. 
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Figures  1 and  2 indicate  much  stronger  near  source  coupling  (in  the 
saturated  clay  and  high  sonic  velocity  layers)  from  the  Baneberry  event  than  from 
a standard  Yucca  Flat  shot.  This  result  agrees  with  previous  analysis  of  Baneberry. 
However,  Baneberry  did  not  couple  as  strongly  as  events  on  Rainier  Mesa  which 
were  completely  contained. 

Near  source  data,  such  as  in  Figs.  1 and  2 give  little  insight  into 
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coupling  at  teleseismic  distances.  J.  Savino  (1975)  of  S has  analyzed  teleseismic 
data  from  Rainier  Mesa  and  Yucca  Flat  and  concludes  that  body  wave  magnitudes 
from  Baneberry  (a  shot  located  above  the  water  table)  were  a factor  of  5 greater 
than  for  other  Yucca  Flat  shots  above  the  water  table,  and  a factor  of  2 greater 
than  Yucca  Flat  shots  below  the  water  table.  This  must  be  due  to  the  saturated 
clay  below  and  at  the  working  point.  However,  it  should  be  noted  that  Rainier  Mesa 
shots  give  body  wave  magnitudes  comparable  to  BancLerry. 

It  is  proposed  here  that  Baneberry  vented  due  to  the  absence  of  an  ef- 
fective compressive  residual  stress  field.  Significant  residual  stresses  are  ob- 
served by  Rimer  (1976)  for  most  of  the  calculations  reported.  Figure  3 shows  the 
compressive  residual  stress  field  around  the  cavity  for  the  Rainier  Mesa  calcula- 
tion discussed  above.  An  overburden  pressure  of  56  bars  must  be  added  to  the 
stress  shown.  These  stresses  tend  to  resist  the  cavity  pressure  and  thus  contain 
the  cavity  gases.  Duff,  et  al.  (1975),  have  reported  that  for  a residual  stress  field 
comparable  to  the  Rainier  Mesa  calculations,  cavity  pressures  almost  three  times 
overburden  can  be  developed  before  any  possibility  of  teneUe  failure  is  to  be  ex- 
pected. Thus,  residual  stresses  result  in  a strong  cavity  structure,  acting  to  pre- 
vent tensile  hydrofracturing. 

No  comparable  residual  stress  field  is  seen  in  the  Baneberry  calcula- 
tion. In  fact,  in  order  to  find  any  residual  stresses  of  significant  magnitudes  in 
Baneberry,  one  must  look  at  the  third  (alluvium)  layer.  Figure  4 shows  the  radial 
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Radial  and  tangential  residual  stress  versus  final 
position  for  Rainier  Mesa  material  properties. 
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stress  distribution  vs  time  at  a range  of  1.20S  x 10  cm.  This  range  is  the  loca- 
tion of  the  largest  radial  residual  stress  at  late  times  and  corresponds  to  a peak 
stress  from  Fig.  2 of  215  bars.  Stresses  are  only  plotted  for  late  times.  Fig- 
ure 5 gives  transverse  stress  vs  time  at  this  range  while  Fig.  6 shows  radial  ve- 
locity at  this  range  decaying  to  zero. 

A comparison  of  residual  stresses  for  the  three  calculations  of  in- 
terest here  is  presented  in  Table  2.  Note  that  the  Rainier  Mesa  calculation  used 
a less  realistic  source  description  which  tends  to  underestimate  these  residual 
stresses. 


Table  2.  Peak  Residual  Stresses  (bars)  for  Baneberry,  Rainier  Mesa 

and  Yucca  Flat 


Radial  Stress 


Transverse  Stress 


Baneberry 

Rainier  Mesa  saturated  tuff 
Yucca  Flat  dry  tuff 


The  Baneberry  calculation  showed  far  lower  residual  stresses  in  both 
the  radial  and  transverse  directions.  This  is  not  surprising  since  the  Baneberry 
working  point  layer  is  believed  to  have  negligible  yield  strength  due  to  the  high 
saturated  clay  content  as  evidenced  by  many  caveins  during  drilling.  Rimer  (1976) 
has  shown  the  magnitudes  of  residual  stresses  to  be  proportional  to  yield  strength. 


TRANSVERSE  STRESS  C 1 0 DYNES/CM**2  D 


VELOCITY  C 1 0 


The  negligible  residual  stresses  seen  for  the  Baneberry  calculation 
can  hardly  be  expected  to  contain  the  cavity  gases.  In  fact,  in  light  of  the  results 
for  the  Rainier  Mesa  calculations,  in  the  absence  of  a significant  residual  field  at 
late  times,  tensile  hydrofracturing  would  be  likely.  No  evidence  is  presented  here 
showing  how  far  from  the  source  fracturing  may  extend  or  whether  it  would  indeed 
extend  to  the  surface.  This  is  a complex  question  involving  the  condensible  flow 
of  steam  in  a growing  crack.  However,  the  paleozoic  escarpment  probably  was 
incidental  to  the  venting,  but  it  may  have  influenced  the  direction  and  propagation 
of  the  crack. 


4.  CONCLUSIONS 

The  calculations  presented  here  do,  in  fact,  show,  as  reported  by 
others,  that  Baneberry  does  couple  higher  than  other  shots  in  Yucca  Flat,  due  to 
the  saturated  clay  about  the  working  point  and  the  high  sonic  velocity  above  the 
working  point.  However,  this  coupling  is  lower  than  for  Rainier  Mesa  shots  which 
consistently  are  contained.  Therefore,  it  is  concluded  that  the  higher  coupling  was 
probably  not  the  cause  of  the  containment  failure  of  Baneberry. 

A definitive  two-dimensional  calculation  of  Baneberry  using  the  latest 
equation  of  state  data  has  yet  to  be  done.  This  calculation  should  be  made  and  should 
incorporate  the  most  sophisticated  treatment  of  tension  failure  possible,  such  as  the 

3 

model  now  available  at  S . Such  a calculation  would  address  the  question  of  whether 
the  impedance  mismatch  at  the  high  sonic  velocity  layer  would  greatly  increase  the 
likelihood  of  tension  failure.  It  is  possible  that  the  tension  failure  induced  by  the  im- 
pedance mismatch  could  destroy  the  residual  stress  membrane  around  the  cavity. 

In  the  absence  of  the  two-dimensional  calculation,  the  following  sce- 
nario is  presented  as  leading  to  the  venting  of  cavity  gases  to  the  surface  to  late 


times.  The  working  point  material  was  extremely  noncompetent  due  to  the  satur- 
ated clay  content,  leading  to  numerous  caveins  during  excavation.  Therefore,  this 
material  would  be  expected  to  have  a negligible  yield  strength.  Since  the  magnitudes 
of  the  compressive  residual  stresses  around  the  cavity  are  proportional  to  the  yield 
strength,  these  stresses  would  be  negligible  for  Baneberry.  In  fact,  the  Baneberry 
calculation  showed  them  to  be  an  order  of  magnitude  lower  than  for  a comparable 
Yucca  event.  Thus,  no  "mystical  magical  membrane"  would  be  developed  around 
be  cavity.  In  the  absence  of  this  membrane  tensile  hydrofracturing  would  be  likely 
to  occur.  If  this  fracturing  was  sufficiently  great,  venting  would  occur.  The  pres- 
ence of  the  paleozoic  escarpment  is  believed  incidental  to  the  venting,  though  it  might 
have  some  influence  on  the  location  of  the  fissure. 

In  conclusion,  the  negligible  yield  strength  of  the  working  point  ma- 
terial was  the  most  significant  factor  in  the  containment  failure  of  Baneberry.  High 
water  content,  while  enhancing  ground  shock  coupling,  was  not  a major  factor.  To 
insure  a structurally  strong  cavity  capable  of  containing  the  cavity  gases,  the  work- 
ing point  should  be  located  in  a competent  material  with  significant  yield  strength. 


REFERENCES 


1.  Bache,  T.  C. , J.  M.  Savino,  N.  Rimer  and  W.  K.  Savage,  "A  Preliminary 
Comparison  of  the  Teleseismic  Magnitude  of  Explosions  at  Rainier  Mesa  to 
Those  in  Dry  Tuff  at  Yucca  Flat,  " S3  Report  SSS-R- 75-251 7 (submitted  to 
DNA),  December  1974. 

2.  Butkovitch,  T.  R. , "A  Technique  for  Generating  Pressure  Volume  Relation- 
ships and  Failure  Envelopes  for  Rocks','  L3  Report  UCRL-51441,  1973. 

3.  Cherry,  J.  T. , J.  Sweet,  and  E.  J.  Halda,  "The  Relation  between  Fracture 
Induced  Porosity,  Site  Geology,  and  the  Containment  of  an  Underground 
Explosion, " S3  Report  SSS-IR- 74-2142,  June  1974. 

4.  Cowperthwaite,  M.,  W.  H.  ?wisler,  "Improvements  and  Modifications  to  the 
Tiger  Code,"  SRI  Report,  January  1971. 

5.  Crowley,  B.  K.,  and  J.  T.  Cherry,  private  communication,  1971. 

6.  Duff,  R.  E. , K.  H.  Lie,  and  W.  7.  Savage,  "Residual  Stress  Fields  in  the 
Vicinity  of  an  Explosion-Produced  Cavity  and  Related  Topics, " S3  Report 
SSS-R-75-2533  (submitted  to  DNA),  June  1975. 

7.  Germain,  L. , "Distribution  of  Physical  Properties  of  Alluvium  and  Tuff  in 
Yucca  Valley,"  L3  Memorandum,  October  23,  1974. 

8.  Hastings,  C. , private  communication,  1975. 

9.  Laird,  D. , private  communication,  1975. 

10.  Ramspott,  L. , private  communication,  1975. 

11.  Rimer,  N. , J.  M.  Walsh,  and  J.  T.  Cherry,  "Teleseismic  Body  Wave 
Magnitudes  as  Predicted  from  Near  Field  Ground  Motion  Calculations  for 
Four  Nevada  Events, " S3  Report  SSS-R-75-2464  (submitted  to  DNA),  October 
1974. 

12.  Rimer,  N. , "The  Relationship  Between  Material  Properties.  Residual 
Stresses,  and  Cavity  Radius  Due  to  a Nuclear  Explosion."  a Topical 
Report,  SSS-R- 76-2907,  submitted  to  DNA  May  1976. 


19 


REFERENCES  (Cont'd) 


13.  Riney,  T.  D. , J.  K.  Dienes,  G.  A.  Frazier,  S.  K.  Garg,  J.  W.  Kirsch, 
D.  H.  Brownell  and  A.  J.  Good,  "Ground  Motion  Models  and  Computer 
Techniques, " S3  Final  Report  3SR-1071  (submitted  to  DNA),  April  1972. 

14.  Savino,  J. , private  communication,  1975. 

15.  Stephens,  D.  R. , H.  C.  Heard,  and  R.  N.  Schock,  "Preliminary  Equation 
of  State  Data  for  Baneberry  Altered  Tuff, " L3  Memorandum  UCID-15874, 
July  13,  1971. 

16.  Tillotson,  J.  H. , "Metallic  Equations  of  State  for  Hypervelocity  Impact," 
General  Atomic  Report  GA-3216,  July  1962. 


20 


* • * r%J  » - 


APPENDIX 


The  Baneberry  calculation  models  the  rock  using  an  equation  of  state 

g 

of  the  form  developed  by  Tillotson  (1962)  and  fit  at  S to  available  shock  data  for 

a number  of  geological  materials.  For  compressed  states  (p<  pQ)  and  for  cold 

expanded  states  (p  < p and  e < e ).  the  pressure  is  given  by 

U 8 

2 

ep  + Ap  + Bp 


For  expanded  states  (o  < o ) where  e > e \ the  pressure  is  given  by 

u s 


The  phase  transition  from  liquid  to  vapor  (p<  p 
mated  as 


0 


and  e < e < e ’)  is  approxi- 


f 


(e-e)P  + (e  ' - e ) P 
' s v ' s s 


e = specific  internal  energy  as  the  material  is  brought  to  vaporization 
s 

temperature 

e ' = additional  energy  required  to  change  the  material  from  a liquid  to 

o 

a vapor  state 

3 

Air- filled  porosity  * is  included  through  the  S porous  crushup 
model  (p-v  model)  in  which  the  pressure  of  the  porous  material  is  described  by 


where 

v = specific  volume  of  porous  material 
e = specific  internal  energy  of  porous  material 
p = pressure  obtained  fron^fequation  of  state 
« = distension  ratio  defined  by  a = v/  v > 1 

where  v is  the  specific  volume  of  the  material  with  zero  air-filled  voids. 

The  distension  ratio  is  required  to  decrease  from  an  initial  value  (at 
p = 0)  down  to  1.0  as  the  pressure  increases  to  pc>  the  pressure  limit  at  which  all 
air-filled  porosity  is  irreversibly  removed.  The  pressure  limit  for  the  reversible 
portion  of  the  pore  collapse,  p , locates  the  boundary  between  two  functions  a (v), 
elastic  and  plastic,  which  together  define  the  crush  curve. 

The  material  strength  model  used  in  SKIPPER  requires  that  the  prin- 
cipal stress  be  within  the  von  Mises  yield  surface.  For  spherical  geometry,  this 
reduces  to  the  condition  that  the  magnitude  of  the  deviatoric  component  of  the  stress 
in  the  radial  direction  not  exceed  2/3  the  yield  strength.  The  yield  function  Y, 
considered  to  depend  both  on  pressure  and  energy,  for  e < e is  given  by 
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- Yo  + YmF-  (2-F-)l  (>-r-) 

m ' nr  J ' m ' 


Y = (Yn  + Y ) 
o m 


M 

' tn  ' 


P<  P. 


P > Pr 


where 


Yq  = yield  strength  at  zero  pressure 

p = pressure  at  which  the  maximum  stress  difference  under  triaxial 
m 

compression  (Y^  + Y ) is  attained 

e = energy  of  melting 
m 

For  energies  greater  than  e , no  deviatoric  stresses  are  permitted. 

m 

Table  A-l  summarizes  the  equation  of  state  parameters  and  elastic 
properties;  sonic  velocity  (c^),  bulk  modulus  (k^),  and  shear  modulus  (G),  used 
in  the  calculations. 
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